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Abstract 

Recently, the search for Majorana fermions (MFs) has become one of the most important and exciting issues in condensed 
matter physics since such an exotic quasiparticle is expected to potentially give rise to unprecedented quantum phenomena whose 
functional properties will be used to develop future quantum technology. Theoretically, the MFs may reside in various types of topo¬ 
logical superconductor materials that is characterized by the topologically protected gapless surface state which are essentially an 
Andreev bound state. Superconducting doped topological insulators and topological crystalline insulators are promising candidates 
to harbor the MFs. In this review, we discuss recent progress and understanding on the research of MFs based on time-reversal- 
invariant superconducting topological materials to deepen our understanding and have a better outlook on both the search for and 
realization of MFs in these systems. We also discuss some advantages of these bulk systems to realize MFs including remarkable 
superconducting robustness against nonmagnetic impurities. 
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I. Introduction 

Majorana fermions (MFs) are neutrinos whose property is 
that particles are their own anti-particles. The existence of MF 
was initially predicted in high energy physics, but has not been 
experimentally proved yet fl|,|2|] . Recently, it has been predicted 
that the MFs can emerge in topological superconductor (TSC) 
as gapless excitations [3i |4jj that are not electrons or holes but 
Bogoliubov quasiparticles formed by bound state due to the su¬ 
perposition of electrons and holes |2|]. The gapless surface or 
edge state within the bulk superconducting gap, regardless of 
whether or not it is fully gapped or partially gapped with nodes 

II , is essentially an Andreev bound state (ABS). It is topologi¬ 
cally protected since its quantum-mechanical wave function has 
nontrivial topology. After the theoretical prediction, the search 
for the MFs in condensed matter physics has become of g reat 
interest 
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2 l 313] because, if discovered, MFs 
will deepen our understanding of quantum states of matter in 
physics and foster innovations in future quantum technologies. 
Furthermore, the discovery would provide us a clue to answer 
what made it difficult for us to observe MF as a neutrino. 

Since recent progress on research of MFs elucidates that the 
topological superconducting states which harbor MFs can be 
realized in many different systems [2j,|<S, ;ll Si i§| E 11, 12, 
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it would be useful to look over the prerequisites for realiz¬ 
ing TSCs. In this review, we focus on bulk time-reversal- 
invariant superconductors (SCs) derived from degenerately- 
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doped topological materials that are considered to be a promis¬ 
ing platform, namely superconducting doped topological insu¬ 
lators (STIs), superconducting doped topological crystalline in¬ 
sulators (STCIs), or superconducting heterostructure materials 
based on doped topological insulators (TIs). We will describe 
advantages/disadvantages of those systems that can be distin¬ 
guished from others in the properties of normal states or sym¬ 
metry invariance. In particular, it is important to notice that 
theoretically the superconducting state in the superconducting 
doped topological materials is robust against non-magnetic im¬ 
purities [28, E T(I 311. Hence, the doped non-magnetic ele¬ 
ment in topological materials will exert almost no ill effects on 
the realization of a TSC with odd-parity paring state. This is 
one of advantages of superconducting doped topological mate¬ 
rials making them to be promising candidates for the TSCs. 

We discuss the experimental efforts and achievements 
on the search for MFs in superconducting doped topologi¬ 
cal materials, candidate materials for time-reversal-invariant 
TSCs: Cu v Bi 2 Se 3 , Sni_,In,Te, (Pbo. 5 Sno. 5 )i_ x In v Te, and 

Cu A (PbSe) 5 (Bi 2 Se 3)6 and list the issues to be addressed in fu¬ 
ture studies. 

2. Superconducting topological materials 

TSCs are accompanied by quasiparticles that have nontriv¬ 
ial topological properties defined in the momentum space. The 
quasiparticles are essentially MFs which are responsible for 
non-local correlation, non-abelian statistics, and Ising magnetic 
response SUSS H. One of the concrete examples of TSCs is 
the superfluid 3 He, which has offered a prototypical system to 
realize various kinds of t opologic al phases, owing to the rich¬ 
ness of symmetry 124, 25, 32, 3.3.0 ■ As for SCs, only a few bulk 
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materials are known as strong candidates for the host of MFs, 
e.g., Sr 2 Ru 04 |34| and some heavy fermion SCs, in particular 
UPt 3 135', 31011. 

One of the keys to realize TSCs and MFs in bulk materials is 
odd parity pairing that satisfies 


relevant to TIs: 


PA(k)P [ = -A(-Jfc), 


( 1 ) 


where P is the inversion operator. This indicates that for a 
single-band system, only spin-triplet pairing can be a candi¬ 
date for bulk centrosymmetric TSCs and MFs. A sufficient cri¬ 
terion for realizing TSCs in odd-parity pairing was derived in 
Refs. HH3 ITTl l39ll . where the topological property is deter¬ 


mined by the number of Fermi surfaces that enclose the time- 
reversal invariant momenta. 

Carrier-doped TIs and topological crystalline insulators 
(TCIs) which we mainly consider in this paper can offer a 
promising platform to realize TSCs and MFs for the follow¬ 
ing two reasons: First, owing to the orbital degrees of freedom 
of electrons embedded in TIs, odd-parity pairing that satisfies 
Eq. © can be realized even in an 5-wave channel of Cooper 
pairs, as a spin-triplet 5-wave orbital-singlet, which can be fa¬ 
vored by phonon-mediated pairing mechanism. All the odd- 
parity pairings in doped TIs satisfy a sufficient criterion for re¬ 
alizing TSCs ll9i[I(3l. Tlll[39ll . Second, contrary to conventional 
wisdom, it has been predicted that a strong spin-orbit coupling 
makes bulk odd parity superconductivity robust against non¬ 
magnetic disorder 1 28 . 1ME1. Hence, carrier-doped TIs can 
expand the horizons of the research field of topological super¬ 
conducting materials. 

In this section, we give an overview of TSCs and time- 
reversal-invariant TIs as their host materials. Special focus is 
placed on Cu v Bi 2 Se 3 as a promising candidate for TSCs. We 
also make a brief overview on the robustness of bulk odd-parity 
pairing against disorders |28, 3(1 31]. In this paper, we intro¬ 
duce the Pauli matrices in the spin, orbital, and particle-hole 
spaces, 5 /y , <x ;( , and r ;( , respectively. The repeated Greek indices 
imply the sum over x, y, z and we set h = 1. 

2.1. ’Ll TIs and Dirac fermions 

To clarify the topological properties of the parent material 
of STIs, we here start to summarize the discrete symmetries 
relevant to TIs. Since many time-reversal-invariant TIs hold 
the inversion symmetry as well as the time-reversal symme¬ 
try 0, the Bloch Hamiltonian, r Hn(k), satisfies the following 
relations. 


WrtW^ = 'Hw(-k), 
TK Tl (k)T~ l = 'Wn(-k), 


( 2 ) 

(3) 


where P = cr x and T = is y K are the inversion and time-reversal 
operators, respectively (K is the complex conjugation opera¬ 
tor). The symmetry ([3]) indicates that the eigenstates of the 
Bloch Hamiltonian, \iinfk)), have Kramers degenerate part¬ 
ners, \u<?\-k)) = T\u^\k)). 

It is further natural to suppose that 'Hjfk) holds a mirror 
symmetry and /V-fold rotation symmetry about the z-axis, be¬ 
cause the discrete symmetries originate in crystalline symmetry 


M'PI\[(k)M = ( H\\(-k x ,k y ,k z ), 
UNTT-nikdlJ^ - 


(4) 

(5) 


Without loss of generality, the mirror reflection plane is set to be 
normal to the T-axis, where the mirror operator M = is x changes 
k and s to (—k x ,k y ,k z ) and (s x ,-s y ,-s z ). The SU(2) matrix 
Un — exp{-icps z /2 ) describes the A-fold spin rotation about the 
z-axis by an angle ip = 27 i/N (N eZ) and Ry is the corresponding 
SO(3) rotation matrix. 

Owing to Kramers degenerate band structure, a minimal 
model for time-reversal-invariant TIs is a 4 x 4 hermitian ma¬ 
trix. A generic form which holds the symmetries ([2} and © can 
be expanded in terms of the five Dirac y-matrices as 'Hjiik) = 
do(k) + Yjjdj{k)yj £1,42]. Following Ref. £cj], we take the y- 
matrices as (y 1 , y 3 , y 3 , y 4 , y$) = (cr x , cr y , o~ z s x , cr z s y , cr z s z ). The 
symmetries © and © guarantee that all the coefficients are 
real and that d a (k) is even on k for a - 0,1 and odd otherwise. 
The discrete symmetries © and © further constraints 'Pl\ \(k) 
in the lowest order on k as 


'Hti (k) = c(k) + m(k)cr x + v z k z cr y + v(k x s) z cr z . 


( 6 ) 


where c{k) - cq + c\k 2 + + k 2 ) and m(k ) = m$ + m\k 2 

ni 2 (k 2 +k 2 ). The Hamiltonian © describes the low-energy band 
structure of topological materials, including BbSe? (TI) |]43 


44], SnTe (TCI) 145. 46[ 47]. We also notice that as the lowest 
order correction to Eq. ©, the effect of the “war ping ” term, 
(k\ + k 3 _)cr z s z , was discussed in TIs [48] and SCs |49], where 
k± = k x + ik y . 

The topology of the parent material is characterized by the 
Z 2 invariant, v 0 0, which originates in the Berry phase 


of \u„ (k)). The inversion symmetry © simplifies the formula 
for the Z 2 number as (-l) y = fj; where (-TA,) = +1 is an 
eigenstate of P at time-reversal-invariant momenta A, that sat¬ 
isfy [7~, 7Tti(A,)] =0. For the form of the Hamiltonian ©, the 
topological invariant v is nontrivial (odd) when sgn(mo'»i) < 0. 
Even in the case that v is trivial, however, another topological 
number, the Chern number, can be introduced in a mirror sym¬ 
metric plane of the Brillouin zone ll5ll 0]. A nonzero mirror 
Chern number defines topological crystalline insulators. As 
will be discussed in Sec. 13.21 the concrete example is SnTe, 
which is accompanied by gapless states on the (001) surface, 
even though the Z 3 number is trivial. 

One of the remarkable consequences of an odd v is the 
emergence of topologically protected Dirac fermions that are 
bound at the surfaces. The wave function of Dirac fermions 
on the surface (z = 0) is obtained by solving the equation, 
[ < H-n(k x ,ky, -id z ) - p]<p D (z) = s(k x ,k y )<p D (z), with a boundary 
condition (fi\}(0) = 0. The wave function with the linear disper¬ 
sion, e(k x , k y ) = ±v Jk 2 + k y , is given for c\ -C 2 = 0 as 


Vviz) = {e~ K ~ z - e~ K + z )\ ” ) »u s (k x ,k y ), (7) 


where k± = v z /2m\± J(mo + t^kp/mi + (v z /2m{) 2 and (0,1)J 
is the spinor in the orbital space ( A T denotes the transpose of 
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a matrix A). The spinor u s in the spin space is determined by 
( k x s y - kyS x )u ± = ±k\\u±. The penetration depth of the Dirac 
cone, {, is defined as ( = kZ 1 ■ The surface Dirac fermions in 
Eq. ([7} are fully polarized in the orbital space. In Sec. 14.11 we 
will discuss the effect of the orbitally polarized Dirac fermions 
on bulk superconductivity. 


2.2. Topology in superconducting topological materials 

To derive a sufficient condition for realizing topological su¬ 
perconductivity in doped topological materials, let us here 
overview the topological invariants, Z and Z2 numbers, rele¬ 
vant to time-reversal-invariant SCs. The topological invariant is 
usually introduced by the global structure of the wave functions 
of quasiparticles that have a fully gapped excitation spectrum. 
Such a gapped SC with time-reversal symmetry is cate gor ized 
to the class Dill in the Altland and Zirnbauer (AZ) table I hi 521, 


and the topologically non-trivial property in three-dimension is 
characterized by a Z topological number, i.e., a winding num¬ 
ber. As clarified in Ref. 113911 . however, topological invariant 
can be introduced even in a SC having nodal structures. In this 
case, the Z is an ill-defined number due to an ambiguity for re¬ 
defining the superconducting gap, while a Z 2 that corresponds 
to the parity of the winding number remains as a well-defined 
number. 

The quasiparticle states that characterize topological super¬ 
conductivity are obtained by solving the Bogoliubov-de Gennes 
(BdG) Hamiltonian in the particle-hole space, 


<U(h\ - I ^Ti(k)-p -iA(k)Sy 

’ 1 iSyA'(k) -Tl^(-k) + p 


( 8 ) 


where p is the chemical potential. In superconducting states, the 
particle sector is coherently coupled to the hole sector through 
the pair potential A (k). Here, we concentrate our attention on 
time-reversal-invariant SCs that satisfy T A(k) r J 1 = A(-k). 
The BdG Hamiltonian, then, holds the time-reversal symme¬ 
try, T7T(k) r ) 1 = 'H(-k), and the particle-hole symmetry, 
C‘H(k)C~' — -'H(-k) with C - t x K. Combining these sym¬ 
metries, one can define the chiral operator T, which is anti- 
commutable with T~l{k), 


ir,9d(k)} = 0, T = -iCT. 


(9) 


The discrete symmetries with T 1 = —1, C 2 = +1, and T 2 = +1 
categorize this system to the Dill class in the AZ table nm. 

The topological properties of superfluids and SCs are gen¬ 
erally determined by the global structure of the Hilbert space 
spanned by the eigenvectors of the occupied band, | u„(k)) ob¬ 
tained from Tl(k)\u n (k)) = E n {k)\u n {k)). To capture the topo¬ 
logical property, it is convenient to introduce the so-called Q- 
matrix, Q(k) = Ze„>o \u n (k))(u n (k)\ - £ E „<0 \u n (k)){u n (k)\ B 
15211 . Since the traceless Q-matr ix satisfies the conditions, 
Q 2 (k ) = +1, the Q -matrix continuously flattens the eigenvalues 
of 'H(k) to -1 (occupied) and +1 (unoccupied). Therefore, the 
Q-matrix maps the Brillouin zone onto the target space spanned 
by the eigenvectors of the BdG Hamiltonian. The chiral sym¬ 
metry is still preserved by the Q-ma trix, {T, Q(k )} = 0, which 
is crucial for introducing topological invariant. In the basis that 


T is diagonalized to T = diag(+l,-1), the Q-matnx becomes 
off-diagonal and is reduced to q(k) e\J{N). Hence, the relevant 
homotopy group for the projector Q(k) in three dimensions is 
given by 713 [U(A)] = Z. The topological invariant is defined as 
the three-dimensional winding number [ 


M |5J, 


W 3d 


f 


dk 

48 ^? 


C /J vri 


tr 


[n Qd M Q)(Qd v Q)(Qd„Q)\ ■ (10) 


For the Dill class, W 3 d can be an arbitrary integer value, which 
guarantees the existence of a surface Majorana cone IS b. BI¬ 
AS clarified in Refs. [Toi, TTj, 39]. the Z 2 number that char¬ 
acterizes the parity of W 3 d determines a sufficient criterion for 
realizing time-reversal-invariant TSCs. Here, we suppose an 
odd-parity pairing that satisfies Eq. 0. This ensures that the 
BdG Hamiltonian holds the inversion symmetry up to the n/2 
phase rotation. 


PTilk)^ = Ti(-k), 


( 11 ) 


where V-Pt : . The inversion symmetry simplifies the parity of 

w 3d to ESInlBii 

(_l) l *' 3 ti - n sgn [£ 2m(A ;)], (12) 

i,m 

where s m (A,j is the energy dispersion of the normal state 
at time-reversal invariant momenta, which is obtained as the 
eigenvalues of 7Tn(A/)-p with the chemical potential p. Owing 
to the time-reversal symmetry, the energy band of carrier-doped 
TIs has a Kramers pair as £ 2m (A,) = £ 2m+ i (A,j. 

The Z 2 number in odd parity SCs was first introduced by 
Sato |Q] in the case of a single-band system and extended to 
a general case by Sato mb and Fu and Berg 10], indepen¬ 
dently. Equation (IT2l) indicates that the mod-2 winding number 
is determined by counting the number of Fermi surfaces en¬ 
closing k = A/. The topological invariant is non-trivial, i.e., 
(-l) li ' 3d = -1, when the number of the Fermi surfaces is odd, 
and trivial otherwise, (-l) M ' 3d = +1. Hence, odd-parity pairing 
and the shape of the Fermi surfaces at time-reversal-invariant 
momenta are sufficient conditions for realizing topological su¬ 
perconductivity and helical MFs in doped TIs. 

Both the Z and Z 2 numbers require the quasiparticle excita¬ 
tion to be fully gapped, i.e., dct'H(k) + 0. For time-reversal- 
invariant nodal SCs, however, the node structure is not topo¬ 
logically protected and removed by introducing a small per¬ 
turbation 5 that transforms A to open a finite gap at the nodal 
points without breaking the time-reversal symmetry J30. The 
unambiguously-defined 6 is set to be 6 —> 0 after the Z 2 (not 
Z) topological invariant is evaluated. This simple procedure 
enables to define the topological numbers in nodal SCs. 


2.3. Topological odd parity superconductivity in doped TIs 
(TCIs) 

Since orbital degrees of freedom is inherent to the parent 
material of STIs, a Cooper pair potential realized in STIs has 
orbital degrees of freedom in addition to the spin and mo¬ 
mentum dependence. Therefore, a general form of the pair 
potential in bulk STIs is given by (/A,v y ) V|V2iCr|Cr2 (A:), where 
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(si, .s’ 2 ) and (<xi, a 2 ) refere to spin and orbital indices of paired 
particles, respectively. In general, the Fermi-Dirac statis¬ 
tics requires the pair potential to be odd under the exchange 
of relative coordinates, and orbital states, (iAs y ) SlS2t0 - l0 - 2 (k) = 
—(iAsy)s 2Sua - 2(ri (—k). For s-wave channel, this condition re¬ 
duces to 


A = j v A s v . 


(13) 


Here, we suppose even-frequency superconductivity as the bulk 
pairing of STIs fl]. In single-band SCs, the condition is satis¬ 
fied only by spin-singlet even parity pairing states. Owing to 
orbital degrees of freedom inherent to the parent material of 
STIs, however, odd parity pairing can be realized even in s- 
wave channel, when the conditions 0} and (IT3l > are satisfied. 

Let us now specify the possible ,v-wave pairing symmetry 
that satisfies the conditions CD and (flTl i for realizing odd par¬ 
ity superconductivity. For spin singlet pairing, the pair poten¬ 
tial A that satisfies Eq. (IT3l> is generally parameterized with a 
vector in the orbital space as A = icr^CTyd^, which is the 
spin-singlet even-parity orbital-triplet state. Among the possi- 
ble orbital-triplet states, only the pairing, A = kr x a y d x oc cr z , 
can be odd parity pairing. A general form of spin triplet s- 
wave pairing is obtained as A = iSyS^d®, with a vector in the 
spin space, d (t \ This pairing symmetry categorized to the spin- 
triplet even-parity orbital-singlet state always satisfies the odd- 
parity condition in Eq. ([[}. 

In Table Q] we summarize possible ,v-wavc pairing symme¬ 
tries in STIs with the rhombohedral D\,i crystalline symme¬ 
try I IS 39, which are relevant to both Cu A Bi 2 Se 3 and 

Sni_ A In A Te (See Sec. 13.11 and 13.2b . There are six indepen¬ 
dent matrices, (A[ a , A\b<r x , A 2 cr v s z , At,ct z , A^ x cr y s x , A^CySy), 
that satisfy the criterion (IT3l) for the Fermi-Dirac statistics, 
where A ; - (j: la, lb, 2, 3, Ax, and 4 y) denote the amplitude of 
each pair potential. Their gap structures are also summarized 
in Table|T] where (A|„, A|/,<x v ) and A 2 cr y s y are fully gapped and 
the others are accompanied by point nodes 056115711 . We notice 
that full gap can be realized in the E u state (see Sec. 14. lb . As 
mentioned above, all odd parity pairing states can be topologi¬ 
cally nontrivial, regardless of nodal structures, when the time- 
reversal-invariant momentum A = 0 is enclosed by the Fermi 
surface. 

A direct consequence of topological odd parity superconduc¬ 
tivity is the existence of gapless surface states. The surface 
structure and tunnelingspectroscopy of STIs have been studied 
theoretically [55, 58| 59, 60j, 6j]]. The A 2u and E u states in Ta¬ 
ble!]] which are topological odd parity pairing with point nodes, 
are accompanied by zero energy flat bands connecting two point 
nodes that are bound to the surface. We will illustrate in Sec. 14.21 
that the zero energy flat band is protected by the hidden Z 2 sym¬ 
metry that is obtained by combining the mirror symmetry and 
time-reversal symmetry. Based on well-established wisdom in 


unconventional SCs B62L [63], the symmetry-protected surface 
Fermi arc is responsible for a pronounced zero-bias conduc¬ 
tance peak in tunneling spectroscopy. 

The A\ u state is similar to the Balian-Werthamer (BW) state 
realized in the superfluid 3 He which is a concrete example of 
time-reversal-invariant TSCs with a full gap |i!,[T, 24, 25, 321. 
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A h 

A 2 <T y s z 
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- 

- 

Z 

A \ u 
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+ 

Z2 

A -2 u 


point 

- 

+ 

Z2 

E u 

A4yCTy Sy 

full 

- 

- 

z 

E u 


Table 1: Pairing potentials in STIs, gap structures, and their parity under inver¬ 
sion P and mirror reflection M = is x : T denotes the representation of the D 34 
symmetry. Note that a full gap behavior in the E u state is discussed in Sec. rm 
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Figure 1: Energy spectra of the A\ u state: The “cone” shape of the surface 
bound states (blue curve) for (mimo/vJ,m 2 »io/v^) = (-0.59,-0.20) (a) and 
the “caldera” shape for (—0.17, —0.20) (b), where Abuik and C denote the bulk 
excitation gap and Fermi momentum, respectively. 


The BW state is accompanied by the gapless surface state hav¬ 
ing a linear dispersion, which is called the Majorana cone. The 
resultant tunneling conductance of the BW state always shows 
a double-peak structure, rather than a zero-bias conductance 
peak [64 ]. However, it has recently been pointed out 155j, 59] 
that the dispersion of gapless surface states twists at finite mo¬ 
menta k\\ = yjk 2 x + k~ which is parallel to the surface. The ve¬ 
locity of the Majorana cone near k\\ = 0, v sur f, is given with the 


parameters in TTri (k) as [55 


Psurf 


yi +Arhi +Am 2 l fx 1 vA 


2ni\jj~ 


m 0 


(14) 


where fh\ = mom\/v At the critical velocity that satisfies v sur f = 
0 , the topological surface state undergoes a structural transition 
from a “cone” (Fig.Q]a)) to “caldera” shape (Fig. [Tib)) l55ll . 

In Fig.[T[b), additional crossings of the dispersion of the sur¬ 
face state appear around k\\ ~k y . The crossing is attributed to the 
cooperative effect of a remnant of the surface Dirac fermions in 
the parent TI and the surface MFs [55|, 59j|- Let us suppose the 
situation that the Dirac cone is well isolated from the conduc¬ 
tion band, which corresponds to the situation with low /.(values. 
This is relevant to CuvBLSej, where well-defined Dirac cone 
was observed by angle resolved photoemission spectroscopy 
(ARPES) [65] at the Fermi level. Therefore, the low energy 
structure is regarded as an effective two-channel model that is 
composed of the conduction band electrons and Dirac fermions. 
The conduction band opens a superconducting gap when adding 
A 2 cr y s z , while the dispersion of the surface Dirac fermions re- 
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main gapless near the Fermi energy [58,166]. Hence, the dis¬ 
persion (IT4l) of the MFs emergent around k — 0 must smoothly 
connect to the dispersion of the gapless Dirac fermions. The 
cooperative effect of surface Majorana and Dirac fermions nec¬ 
essarily gives rise to the twisting of surface Majorana cone and 
results in a caldera structure. For a large p where the surface 
Dirac cone merges into the bulk conduction band, the surface 
states induced by topological odd parity superconductivity re¬ 
mains as a cone shape. The structural transition of the surface 
bound states occurs in the case of the E u state as well. 

Owing to the structural transition, a pronounced zero-bias 
conductance peak can be expected not only in the E u state but 
also in the A i„ state. The A i„ state with the caldera shape of 
surface MFs has a sharp zero energy peak in the surface den¬ 
sity of state, which is responsible for the zero-bias conductance 
peak j55 MSS. 


2.4. Effect of disorders on bulk odd parity superconductivity 

As clarified in Sec. 12.31 odd parity pairing is a sufficient con¬ 
dition for realizing topological superconductivity in doped TIs. 
Similarly with a conventional s-wave pairing, the ,4 1 g state in 
Table Q]is highly insensitive to disorder [ 68 j. Contrary to a con¬ 
ventional s-wave pairing, however, unconventional supercon¬ 
ducting states have been believed to be fragile against disorders. 
This suggests that unconventional superconductivity can be re¬ 
alized only in clean materials. Since the superconducting doped 
topological materials typically have a short mean-free path [[690, 
the robustness of odd-parity superconductivity remains an im¬ 
portant issue. 

Michaeli and Fu |[28ll recently predicted that in contrast to 
conventional wisdom, the A i„ state, which is one of the odd 
parity pairings, remains robust against disorder. They found 
that an approximate chiral symmetry emergent in STIs prevents 
the superconducting compound from pair decoherence induced 
by impurity scattering. Similar conclusions were derived in 
Ref. 13011 . based on a self-consistent T -matrix approach for im¬ 
purity scattering. The robustness against disorder was also re¬ 
ported by Nagai 113 ill in the case of the E u state which has the 
odd-parity pairing with point nodes. 


3. Candidate materials 

The search for topological phases of matter with time rever¬ 
sal symmetry brought about the discovery of TI materials and 
stimulated the search for an even more exotic state of matter, 
TSC, a superconducting cousin of TI. Opposed to the ideal TIs, 
the bulk of existing materials at the early stage of TI research 
was conductive and smeared the novel properties of the sur¬ 
face states. While achieving a bulk-insulating state in real TI 
samples has been a great experimental challenge, inducing su¬ 
perconductivity in TI samples by doping carrier proved to be a 
great experimental challenge as well. 

3.1. Electron-doped TI: CuxBTSe 3 

Cu A Bi 2 Se 3 , the first example of a superconducting doped TI, 
is a layered material consisting of stacked Se-Bi-Se-Bi-Se quin¬ 
tuple layers that are only weakly bonded to each other by van 


der Waals forces. It has the same rhombohedral crystal struc¬ 
ture as Bi 2 Se 3 (R3m space group Mlzi) with the representa¬ 
tion of point group [10]]. Its superconductivity arises as a 
consequence of Cu intercalation into the van der Waals gap of 
the parent compound Bi 2 Se 3 yielding an electron concentration 
of ~2xl0 20 cnT 3 with a maximum transition temperature T c of 
3.8 K for 0.12 < x < 0.15, as well as a metallic behavior of 
the resistivity in the normal state s. The bulk superconduc¬ 
tivity of Cu A Bi 2 Se 3 samples was confirmed by the specific heat 
measurements which showed that the superconducting state of 
Cu A Bi 2 Se 3 for x - 0.3 is likely to be fully gapped [[69]. The 
carrier density of order of 10 20 cm -3 is very low for SCs whose 
maximum T c is in the order of a few K within the framework of 
the Bardeen-Cooper-Schrieffer(BCS) theory. The anomalously 
high T c for a relatively low carrier density may indicate that the 
pairing mechanism of the superconducting Cu A Bi 2 Se 3 could be 
affected by the strong spin-orbit coupling, albeit still be gov¬ 
erned by the electron-phonon mechanism in the same way as 
most of the existing SCs. 

Stimulated by the discovery of superconductivity in doped TI 
Cu A Bi 2 Se 3 , in 2010, Fu and Berg discussed the possible pairing 
symmetry within a two-orbital model and a sufficient criterion 
for realizing time-reversal-invariant TSCs in materials with in¬ 
version symmetry Gi. Intriguingly, a novel spin-triplet pairing 
with odd parity is favored by strong spin-orbit coupling with 
gapless surface ABSs within a full pairing gap in the bulk: We 
emphasized in Sec. 2.2 that what is important for TSCs is odd 
parity pairing (See Eq. 0}) and a nontrivial Z 2 number is de¬ 
fined for all odd parity pairing even in nodal SCs. 

The discovery of the Cu A Bi 2 Se 3 SC attracted a lot of atten¬ 
tion in the community, however, experimental studies of this 
material were hindered at the beginning by difficulties of ma¬ 
terial synthesis. Thanks to employing electrochemical Cu in¬ 
tercalation (EC1) technique combined with the post annealing 
that is essential for the activation of the superconductivity (a 
typical sample image of EC1 Cu A Bi 2 Se 3 is displayed in Fig. 
2(a)) 069l 17ill , the maximum shielding fraction of Cu A Bi 2 Se 3 
now can reach up to ~50% at 1.8 K with Cu concentration x of 
~0.3; besides, the maximum value can be even larger at lower 
temperatures. Meanwhile, serious improvements of the melt- 
grown (MG) synthesis have been done 13 3 74j|. Currently 
Cu A Bi 2 Se 3 has been the most widely studied as a candidate for 
the TSCs. 

It would be prudent to discuss disorder effects in Cu A Bi 2 Se 3 
since odd-parity pairing state in general disfavors disorders 
175|, 76, 77] as mentioned in Sec l2.4l There is a systematic 
unbalance between the number of the intercalated Cu atoms 
and that of the increase in the carrier density ll69h . Progres¬ 
sive Cu intercalation into the system introduces si gnif icant dis¬ 
order 170. T\_. 731 and an intrinsic inhomogeneity H7 lh . leading 
to an anomalous suppression of the superfluid density which 
was deduced from the measurements of the lower critical field 
1 78] and an optical spectroscopic study |79[. At the same time, 
the transition temperature T c is monotonously and only mod¬ 
erately suppressed from T c - 3.8 K (x ^ 0.1) to T c - 2.6 K 
{x - 0.5), which agrees with the non-magnetic impurity effect 


theoretically predicted for this class of TSCs 128, [2^, 3(j, 31]. 
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Figure 2: Cu r Bi 2 Se 3 and zereo-bias conductance peak, (a) A typical cleaved 
surface image of Cu-intercalated Bi 2 Se 3 with post annealing, (b) Point-contact 
spectra ( dl/dV vs bias voltage) of Cu v Bi 2 Se 3 with x = 0.3 for 0.35-2 K mea¬ 
sured in 0 T. 


Figure 3: Legitimacy of soft point contact, (a) Atomic force microscope pic¬ 
ture of Ag nano grains on the measured surface after the measurement lead is 
removed, (b) A diameter of the contact area of a nano grain in a Ag cluster at 
the sample surface, (c) An expected diameter dependence in the resistnace of 
an electron-transport channel. 


The robust topological superconductivity even in disorder can 
be derived from the relativistic effects in doped topological ma¬ 
terials I30ljlh. 

In order to confirm whether or not a SC is a TSC, it can be 
useful to verify the existence of an ABS consisting of MFs that 
is a low-energy gapless excitation at the surface by detecting 
surface “low-energy” phenomena, in particular electron tunnel¬ 
ing. The gapless excitations that create the density of states 
within the superconducting gap can be detected as a zero-bias 
conductance peak (ZBCP) in tunneling spectroscopy. Indeed, 
both the A\ u and E u states in Table [T] which have bulk topo¬ 
logical odd parity superconductivity, are responsible for ^pro¬ 


nounced ZBCP on the (111) surface of Cu A Bi 2 Se 3 [55, 59]. It 


should be mentioned that a ZBCP was observed in S^RuO/j, a 
strong candidate for time-reversal breaking chiral p -wave TSCs 


[721 8 1 


Actually, the ZBCP was observed by soft point-contact 
spectroscopy experiments on a cleaved (111) surface of ECI 
CUjBhSej for x - 0.3 with a tiny drop of silver paint con¬ 
sisting of clusters of Ag grains that produce Ag-grain parallel 
channels at the point contact as shown in Fig. 2(b) |3^]. A typi¬ 
cal diameter of grains appears to be ~20-50 nm estimated from 
the image taken by the atomic force microscope in Fig. 3(a) 
and the diameter of the contact area for one grain d is always 
smaller than that of the grain size (see Fig. 3(b)). In the case 
of ballistic electron transport, the resistance of the channel can 
be estimate as Sharvin resistance that is proportional to 1/d 1 
while in the case that electrons are scattered and loose kinetic 
energy, the channel resistance can be estimated as Maxwell re¬ 
sistance that is proportional to 1/d 18111 . The Sharvin resis¬ 
tance becomes dominant when d is smaller than the intersec¬ 
tion point ofj between two types of resistance with respect to d 
as indicated in Fig. 3(c), which is ~ 100 nm in the case of Ag 
182||. Since the Ag grain size is smaller than c/p the electron 
transport in the case of the soft point contact can be ballistic, 
and hence spectroscopy in the presence of a finite repulsive po¬ 
tential barrier at the boundary between the Ag grain and a SC 
can be performed. Theoretical considerations of all possible 
superconducting states shows that the observed surface ABS 
of Cu v Bi 2 Se 3 is due to MFs [39]. Consistently, similar ZBCPs 
were observed by other point-contact spectroscopy experiments 


with Au tip on MG Cu A Bi 2 Se 3 for x - 0.25 

On the contrary, both unconventional ZBCP and conven¬ 
tional Blonder-Tinkham-Klapwijk (BTK) spectra were ob¬ 
served on the surface of MG Cu A Bi 2 Se 3 for x - 0.20 by point- 
contact spectroscopy with Au tip & Andreev reflection spec¬ 
troscopy performed by a nanoscale device showed spectra par¬ 
tially fitted by BTK theory [84] ■ Moreover, a conventional 
tunneling spectrum was observed on a cleaved surface of ECI 
Cu A Bi 2 Se 3 for x - 0.25 by scanning tunneling spectroscopy 
(STS) [85] ■ These aroused a controversy with respect to the na¬ 
ture of superconductivity in Cu A Bi 2 Se 3 . It is worth noting that 
recent self-consistent calculations for surface density of states 
(SDOSs) has clarified this puzzling issue: the surface Dirac 
fermions in a conventional pairing state of the ,v-wave SC will 
open an additional gap larger than the bulk superconducting gap 
yielding a two-gap structure due to parity mixing of the pair 
potential near the surface, nonetheless the STS spectra exhib¬ 
ited only single-gap structure. Importantly, the parity mixing 
is absent in an unconventional pairing state of odd-parity SC 
& In addition, recent ARPES experiments revealed that the 
Fermi surface (FS) of superconducting Cu A Bi 2 Se 3 is a cylindri¬ 
cal open FS [86] coexisting with the topological gapless sur¬ 
face states of Bi 2 Se 3 well-separated from the bulk conduction 
band 65t 8^1 . Therefore the STS spectra can be originated from 
topological odd-parity superconductivity with a cylindrical FS 
B67I1 . Further details will be discussed in Sec. 14. II 
It is worth mentioning that the pressure variation of the su¬ 
perconducting phase of MG Cu A Bi 2 Se 3 for x - 0.3 was inves¬ 


tigated by applying a hydrostatic pressure up to 2.31 GPa [87] 
and the superconductivity appears to be robust up to 6.3 GPa. 
The sufficiently large mean free path and the variation of B C 2 (T) 
that agrees with the p-wave polar-state model indicate the p- 
wave spin-triplet superconductivity of Cu A Bi 2 Se 3 , although an 
anisotropic spin-singlet state cannot be discarded completely. 
Apparently further studies on Cu A Bi 2 Se 3 with different mea¬ 
surement techniques are required to elucidate the true nature of 
this material. 


3.2. Hole-doped TCI: Sn\- x In x Te and (Pbo.sSno. 5 )\- x In x Te 
Based on a guiding principle suggested from the discovery 
of Cu A Bi 2 Se 3 , low-carrier-density semiconductors whose Fermi 
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Figure 4: Sni_ v In v Te. (a) A typical facet image of Sni_ v In A Te crystal grown by 
the vapor-transport method, (b) Cubic rocksalt crystal structure of Sni_ A In A Te. 
(c) The Fermi surfaces of p-type SnTe are centered around the four equivalent 
L points, which belong to the time-reversal-invariant momenta, in the bulk Bril- 
louin zone of the cubic structure with fee Bravais lattice, (d) Point-contact spec¬ 
tra (Bias-voltage (V) dependence of the differential conductance, [dI/dV](V), 
at various temperatures) of Sni_ A In A Te with x = 0.045 measured in 0 T. 


surfaces are centered around time-reversal-invariant momenta 
in strong spin-orbit coupling systems, in particular doped TCIs 
with mirror symmetry, are other likely candidates for TSCs. In- 
doped SnTe (denoted Sn i_ f In A Te), a first example of STCIs, has 
the same rocksalt crystal structure ( Fm3m space group) with 
the representation of O /, point-group symmetry as pristine SnTe 
(see Fig. 4(b)) at ambient temperature. It reveals a metallic 
behavior of the resistivity in the normal state in single crys¬ 
tals synthesized in two different ways: Melt growth [[88] and 
vapor-solid growth [89, 90]. A typical faceted single crystal is 


shown in Fig. 4(a). It has been known that hole-doped SnTe, 
a degenerately-doped TCI due to Sn vacancies, superconducts 
and T c increases as the hole density increases from 0.024 K 
(~5xl0 20 cm -3 ) to ~0.2 K (20x 10 20 cm -3 ) [HI. When cooled, 
ferroelectric structural phase transition (SPT) occurs in SnTe 
so that superconducting degenerately-doped SnTe must have a 
rhombohedral crystal structure. 

Upon substitution of Sn by In, however, T c jumps as high 
as ~1 K and gradually increases as the In content increases, 
suppressing the SPT, although the hole density remains simi¬ 
lar [88]. Above x ~ 0.04 the SPT is completely suppressed. 
Nevertheless, the topological surface state of SnTe is preserved 
on Sni_Jn v Te for x = 0.045 providing evidence for the pres¬ 
ence of an inverted band structure as in pristine SnTe [92] and, 
therefore, Sni_ v In v Te with a moderate In content is a super¬ 
conducting doped TCI. Although it has not been clarified yet 
whether or not Sn _Jn v Te with very high In content can still be 
a doped TCI, interestingly, progressive In doping enhances T c 
up to 4.5 K for x = 0.45 Hj or -4.8 K for x = 0.45 0 de¬ 
termined by the onset of diamagnetism, and 4.7 K for x - 0.40 
determined by the onset of superconducting transition in the re¬ 


sistivity measurement | 

The bulk superconductivity of Sni_ x In x Te was confirmed 
by the specific heat measurement for 0.021 < x < 0.12 [88], 
0.025 <x<~0.05 l%|] and x = 0.4 [95], as well as muon-spin 
rotation or relaxation (//SR) measurement for 0.38 < x < 0.45 
194]]. Both measurements show the superconducting state of 
Sni_ x In v Te is likely to be fully gapped. The volume fraction of 
Sni_ x In v Te is essentially 100%. The absence of impurity phase 
is an advantage of this material over Cu A Bi 2 Se 3 , if this is really 
a TSC. 

It turns out that the superconducting state of Sn i_ A In A Te has a 
rich phase diagram where the dependence of T c on indium con¬ 
tent x below -0.04 (in the ferroelectric rhombohedral phase) 
in vapor-grown (VG) single crystals 1961 appears to be differ¬ 
ent from that in MG single crystals [88]]. The different trend of 
T c can be explained by the unusual role of impurity scattering 
within the framework of the BCS theory based on the Martin- 
Phillips theory for the enhancement of T, [96] implying that the 
strength of the impurity scattering in crystals is probably differ¬ 
ent when grown in different ways. Hence, the superconducting 
state of Sn i_ A In f Te in the ferroelectric rhombohedral phase is 
not necessarily unconventional and can be topologically trivial 
with conventional even-parity pairing. 

Importantly, the ZBCP, signature of the surface ABS 
was observed by the soft point-contact spectroscopy on VG 
Sn i_Jn ,Te single crystals for x = 0.045 as plotted in Fig. 
4(d)[89|, indicating the realization of a topological supercon¬ 
ducting state in the cubic phase. To identify the nature of the su¬ 
perconducting state in Sn i-Jn/fe, the conduction and valence 
bands in the vicinity of four L points of the fee Brillouin zone, 
where ellipsoidal Fermi surfaces in the normal state are located 
as drawn in Fig. 4(c), are described by the k ■ p Hamiltonian 
in Eq. (0 [46]. The effective Hamiltonian of Sit| Jn/I’e has 
essentially the same form and the same symmetry classification 
of the possible gap functions as that of Cu v Bi 2 Se 3 except for the 
presence of the threefold rotation symmetry around (111) axis 
in Sn | Jn v Te. Theoretical considerations lead to the conclu¬ 
sion that all possible superconducting states are topologically 
nontrivial [89]. 

However, even in the cubic phase, no clear bulk unconven¬ 
tionality was observed in Sni_ v In v Te crystals with a very high 
In content so far [93, 94[ 3- It is worth noting that the re¬ 
cent numerical calculations using a self -consistent '/’-matrix 
approach in the case of /.'-independent pairing reveal that the 
superconducting state can be altered from p- (odd-parity pair¬ 
ing) to ,v-wave (even -parity pairing) character depending on the 
magnitude of the relativistic effects in the normal-state Dirac 
Hamiltonian of three-dimensional TSCs [ 3oL 3Tll . In this re¬ 
gard, s-wave superconductivity which is robust against a strong 
impurity scattering would dominate the superconducting state 
of high-In-content samples where p is separated from the Dirac 
points. Hence, the odd-parity state can be realized only in the 
lowest T c samples in the cubic phase where p is close to the 
Dirac points and the impurity scattering is the weakest, though 
the surface states are significantly broadened due to strong 
quasiparticle scattering caused by In doping 192 1. 

Recently the effect of the In concentration on the crystal 
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Figure 5: Specific heat of Cu A (PbSe) 5 (Bi 2 Se 3 ) 6 . Superconducting transition in 
c e i/T in 0 T obtained after subtracting the phonon contribution determined in 2 
T, where c e i is the electronic specific heat. The dashed line is the weak-coupling 
BCS behavior (coupling constant a = 1.76) for T c of 2.85 K. The green solid 
line is the theoretical curve for d -wave pairing on a simple cylindrical Fermi 
surface with line nodes along the axial direction f98h . Horizontal solid line 
corresponds to the normal-state electronic specific-heat coefficient . 


structure and superconducting properties of another supercon¬ 
ducting doped TCI (Pbo. 5 Sno. 5 )i-*In v Te has been investigated 
0. The single crystal of (Pbo. 5 Sno. 5 )i_ A In A Te grown by a 
modified floating-zone method retains the rocksalt structure up 
to the solubility limit of In (x — 0.30). The dependence of T c 
and the upper critical magnetic field (H c 2 ) on the In content x 
has been measured and the maximum T c = 4.7 K for x = 0.30 
with noH c 2 (T = 0) ~5 T was found. Further studies on the ma¬ 
terial to find a signature of unconventional superconductivity 
are necessary to conclude that (Pbo. 5 Sno. 5 )i- A In A Te is a TSC. It 
is also important to search for other candidates for TSCs based 
on superconducting doped TCI. 


3.3. Electron-doped natural heterostructure material: 
Cu x (PbSe) 5 (Bi 2 Se 3 ) 6 


Recently, a new SC based on a topological insula¬ 
tor heterostructure material, Cu v (PbSe) 5 (Bi 2 Se 3)6 (abbrevi¬ 
ated CPSBS), was discovered [99]. The pristine material 
(PbSe)g(Bi 2 Se 3)6 (PSBS) is a natural heterostructure of a TI 
(BFSeg) and an ordinary insulator (PbSe). It was found that the 
PbSe unit works as a block layer and the topological boundary 
states are encapsulated in each BUSe 3 unit, making the sys¬ 
tem to possess quasi-two-dimensional states of topological ori¬ 
gin throughout the bulk GM- CPSBS is synthesized by inter¬ 
calating Cu into PSBS with post annealing that is essential to 
activate the superconductivity similar to the activation process 
of Cu A Bi 2 Se 3 . It is worth noting that the specific-heat behav¬ 
ior of CPSBS suggests for the first time in a Tl-based SC that 
unconventional superconductivity occurs in the bulk with gap 
nodes as indicated in Fig. 5. The existence of gap nodes in 
a strongly spin-orbit coupled SC gives rise to spin-split ABSs 
that are hallmark of topological superconductivity. Hence this 
new SC emerges as an intriguing candidate for TSCs. 



E/ Abuik E/A bulk E /Abuik 

Figure 6: (a) Surface density of states for the bulk s -wave pairing (A]g) at 
(mi,m 2 ) = (-0.17,-0.20). Surface density of states for the A\ u state with a 
spheroidal Fermi surface (b) and a cylindrical shape (c). The insets in (b) and 
(c) show the Fermi surface, where a and c are the lattice constants. 


4. Outlook 


4.1. Pairing symmetry of Cu x Bi 2 Se 2 

Here, we revisit the pairing symmetry and topological super¬ 
conductivity of Cu A Bi 2 Se 3 . As shown in Sec. 3.1, surface sen¬ 
sitive experiments observed pronounced ZBCPs on the (111) 
surface [ 39. 74[. 830, which are strong signatures of bulk topo¬ 
logical odd parity superconductivity. On the other hand, con¬ 
flicting experimental results of conductance and tunneling spec¬ 
troscopy were recently reported [84. El. 110ill, which led to 
contrary statement that this material has a conventional s-wave 
pairing symmetry. 

First of all, we would like to emphasize that, even if the 
bulk of Cu A Bi 2 Se 3 is a conventional .v-wave pairing, i.e., the 
,4 1 g state (Ai,) in TableQ] the surface structure becomes uncon¬ 
ventional [67]- As clarified in Eq. ©, the Dirac fermion is fully 
polarized in the orbital space, which breaks the inversion sym¬ 
metry. Since Cooper pairs are equally populated to both orbitals 
in the A i g state, the orbital polarization of the Dirac fermion 
strongly suppresses the condensation energy at the surface. As 
demonstrated in Ref. [67D, the mixing of the subdominant com¬ 
ponent Agcr,, which has the odd parity pairing, is indispensable 
for the maximum gain of the condensation energy. When the 
surface Dirac fermions that are remnants of the parent TIs are 
well defined, therefore, the intertwining with bulk superconduc¬ 
tivity gives rise to a large energy gap in the surface Dirac cone, 
A sur f. This is distinguishable from the gap in the bulk conduc¬ 
tion band, Abuik- This results in the double-coherence structure 
of the surface density of states, as shown in Fig. [3 a). Hence, 
the bulk conventional even-parity scenario is inconsistent with 
a simple U-shaped form in the surface density of states reported 
in Ref. [85| and the STS result in Sec. l3.ll is a puzzle. 

A possible scenario for Cu A Bi 2 Se 3 is bulk odd parity su¬ 
perconductivity with a Fermi surface evolution. For bulk odd 
parity pairing, no surface parity mixing is induced by Dirac 
fermions. The surface density of states has a pronounced zero 

Us 


energy peak which is responsible for a ZBCP |55j, 58,59 


shown in Fig.[6lb), when the Fermi surface encloses the I point. 
As discussed in Sec. 12.31 the topological superconductivity and 
the existence of zero energy density of states on the surface are 
sensitive to the shape of the Fermi surface. For a cylindrical 
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Fermi surface which does not enclose the F point, the surface 
state is no longer topologically protected and thus the resultant 
surface density of states on the (111) surface becomes a simple 
U-shaped form (Fig.[6jc)). This feature is commonly applicable 
to the Ai„ and E u states. 

We would like to notice that the E u state is the two- 
dimensional representation and the arbitrary linear combina¬ 
tion, cr y s x cos <p+crySy sin </>, rotates the nodal direction, where f 
is the azimuthal angle and the nodal direction is denoted by the 
unit vector (sin <p , cos f). Since the point nodes in the E u state is 
protected by the mirror reflection symmetry ( fl5l i with respect to 
the y-z plane, the bulk excitation in the E u state becomes gapful 
for <p4 0, ±7t/ 3, and ±2n/3. Fu J49] recently suggested that in 
the presence of the extra “warping” term in Eq. ©, the E„ sce¬ 
nario is able to explain all experimental results including both 
the bulk and surface measurements. The fully gapped E u state 
may be energetically competitive to the A\ u state. 

We also mention that in the presence of a phonon-mediated 
short-range interaction, bulk odd parity pairing is energetically 
competitive to even parity pairing (A lx ,). Brydon et al. demon¬ 
strated that in addition to the phonon-mediated interaction, a 
repulsive electron-electron interaction based on the Coulomb 
pseudopotential is critical to stabilizing the spin-triplet odd- 
parity states Il02ll . The bulk odd-parity pairing state, the A |„ 
or fully gapped E u state, may be stabilized by phonon mecha¬ 
nism in the presence of a weak electron-electron correlation. 

Apart from k-independent pairing, mo mentum-dependent 


pairing was examined in Ref. 1 1103 . 104 1. Recently, first 


principle linear-response calculations predict that a p-wave-like 
state can be favored by a conventional phonon-mediated mech¬ 
anism 1110511 . This is attributed to a singular behavior of the 
electron-phonon interaction at long wavelengths. 

4.2. Smoking-gun for topological superconductivity 

As discussed in Sec. 14.11 a promising candidate for 
Cu v Bi 2 Se 3 is an odd parity pairing, the A or E„ state. Since 
both odd parity pairings have a twisted Majorana cone, they can 
equally explain a pronounced ZBCP in tunneling spectroscopy. 
Thus, it is hard to identify the bulk pairing symmetry via point 
contact measurements. Here, we summarize smoking-gun ex¬ 
periments for identifying bulk pairing symmetry and for ex¬ 
tracting properties inherent to topological superconductivity. 

Spin susceptibility. — First, the Knight shift in nuclear mag¬ 
netic resonance experiments can provide a smoking-gun evi 
dence for the spin state of the_bulk pairing symmetry 1561157 


[Toil. Hashimoto et al. l56l [57]] introduced the d-vector that 


describes the spin-triplet component of A in the band represen¬ 
tation. It turns out that the Knight shift in the A |„ state is dis¬ 
tinguishable from that of the E„ state, when an applied field is 
rotated within the ab plane. The former is characterized by an 
hedge-hog-like d-vector profile in the k-space, similarly with 
that of the BW state in 3 He. This is responsible for an isotropic 
magnetic response. In contrast, the E u state results in the uniax¬ 
ial anisotropy of the Knight shift in the ab plane. Furthermore, 
Nagai predicted that nodal structure in the E u state is detectable 
through angle-resolved heat capacity and thermal conductivity 
measurements 1110711 . In the E u state, the zero-energy density of 


states around a vortex core shows twofold rotational symmetry 
and splits along the nodal direction with increasing energy. 

Thermal Hall conductivity. — One of fingerprint experiments 
for bulk topological superconductivity is the quantized ther¬ 
mal Hall conductivity. The topology of three-dimensional fully 
gapped SCs with time-reversal symmetry is characterized by 
the Z topological number, W 3 d, as defined in Eq. (fTTTb . It has 
been shown that the thermal Hall conductivity K xy is quantized 

7[~ k 2 T 

in Z time-reversal invariant TSCs, = i r ,1 W 3 d, when a small 
gap is induced in the Majo rana cone 1 108 . 1091 1 id 111 


11211 . 


Shimizu and Nomura 1113i1 evaluated the thermal Hall conduc¬ 
tivity in the Ai„ state of Cu A Bi 2 Se 3 . They demonstrated that K xy 
depends on a time-reversal breaking perturbation which induces 
a small gap in the Majorana cone and W 3 d is directly detectable 
through the thermal Hall conductivity when the small gap is 
induced by the complex .v-wave pair potential. 

Topological “crystalline” superconductivity .— Topological 
superconductivity peculiar to nodal SCs, i.e., the E u state, is 
unveiled by introducing the combined Z 2 symmetry which is 
composed of the mirror reflection symmetry and time-reversal 
symmetry <[3j. The superconducting state retains the mirror 
symmetry if the gap function A (k) is even or odd under the 
mirror reflection, MA(k)M T = ±A(-k x , k y , k-). Then, the BdG 
Hamiltonian "Hulk) preserves the mirror reflection symmetry 

M ±f H(k)M ±] = ( H{-k x , k y , k z ), (15) 

when an external field is absent. The mirror reflection operator 
AH is defined as Al* = diag(M, ±M*). 

Once the Zeeman fields, H, are applied, the mirror symme¬ 
try with respect to the plane parallel to the Zeeman filed is lost 
and the time-reversal symmetry is explicitly broken, but a com¬ 
bination of them can be still preserved if H ■ x = 0. This is 
because the combination of the mirror reflection and the time- 
reversal rotates the magnetic field H —> (~H X , H y , H z ). Conse¬ 
quently, the Hamiltonian 4T(k) with H x = 0 holds the following 
Z 2 symmetry, 

TM ±, H(k)M ± ^T l = <H(k x , -ky, -k z ). (16) 

Combining the Z 2 symmetry with the particle-hole symmetry, 
C'H z fi (k)C 1 = - r H* s (-k), we define the chiral symmetry oper¬ 
ator, F1 = CT Al*. Then, it turns out that I | is anti-commutable 
with the effective Hamiltonian jr 1 , '71(0, k y , k z ) j = 0. Therefore, 
similarly to W 3 d in Eq. (ITOl) . the one-dimensional winding num¬ 
ber is defined as w\d(k y ) = ~yyj f dkAFEl ^lfodi 'H(k)\i -n lEil 
114 . 115 ]. From the generalized index theorem in Ref. Il 14ll . 
the non-zero value of w\d is equal to the number of zero energy 
states that are bound to the surface. 

Let us now consider A = A^ x a y s x in the E u state, where the 
point nodes lie in the mirror plane = 0). In this situation, 
one finds wid(k v ) = 2 for |k v | < kp, and otherwise wid(k y ) = 0. 
There appear doubly-degenerate zero energy states along the 
chiral symmetric plane k x — 0 in the E u state. Hence, the mirror- 
symmetry-protected topological invariant ensures the existence 
of surface Fermi arc that connects two point nodes. Such 
a surface Fermi arc can also be realized in the planar state 
of 3 He [24] and the £j„ scenario of the heavy fermion SC 
UPt 3 I35D(see Fig. 7). 
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between superconducting TIs, the Josephson current is highly 
sensitive to the difference in the spin-helicity of the surface MFs 
at the interface, where the positive (negative) helicity character¬ 
izes the existence of Majorana “cone” (“caldera”). The Joseph¬ 
son current between different helicity is strongly suppressed at 
low temperatures, which is contrast to that in the junction with 
same helicity. 


5. Summary 


(a) 


(b) 


Figure 7: (a) A stereographic view of the energy gap of the E u state. The Fermi 
surface and point nodes are projected onto the surface in the k x k y plane, where 
the topologically protected surface Fermi arc connects to two point nodes, (b) 
Crystal structure of Cu v Bi 2 Se 3 viewed from the c-axis and three mirror planes. 


Now following Refs. I124. ll 15H116lll 1711 . one can show that 
multiple Majorana zero modes with chiral symmetry ensures 
the Ising character of the topologically protected zero energy 
states. The topologically protected surface Fermi arc does not 
contribute to the local density operator, p (surf ’(r) = 0. This indi¬ 
cates that the MFs protected by the chiral symmetry can not be 
coupled to the local density fluctuations and thus are very robust 
against non-magnetic impurities. Similarly, the local spin oper¬ 
ator, S, is constructed from the surface MF as S (surf) = (S x , 0,0). 
This implies that only the surface MF has anisotropic magnetic 
response. The surface MF and topological Fermi arc are in¬ 
sensitive to a magnetic field along the chiral symmetric plane, 
but fragile against a field perpendicular to the plane. Since the 
surface Fermi arc is responsible for a ZBCP, the anisotropic 
magnetic response of MFs might be observable in point con¬ 
tact experiments when an applied field is rotated within the ab 
plane. The anisotropic magnetic response can be a manifesta¬ 
tion of TSCs with the combined Zi symmetry, i.e., the E u state 
with point nodes. 

We would like to notice that arbitrary linear combination, 
( T y s x cos 4> + cr y Sy sin (f> , cannot be invariant under the mirror re¬ 
flection symmetry M = is x , when + 0, ±n/ 3, and ±2n/3. As 
mentioned in Sec. 14. H and l49fl . the E u state without mirror sym¬ 
metry becomes fully gapped. In this situation, wia introduced 
above is irrelevant and the argument on Majorana Ising spins is 
not applicable to the E u state where the nodal direction is tilted 
from the mirror invariant plane. 

Josephson coupling .— The existence of the Majorana cone 
twisted by the surface Dirac fermions is peculiar to super¬ 
conducting doped topological materials. It has been pre¬ 
dicted IQ, II18 1 that bulk odd parity pairing and surface heli¬ 
cal MF are detectable through Josephson currents. In a junc¬ 
tion between a conventional s-wave and the A\ u state of the 
superconducting doped topological materials, the leading order 
of the Josephson current J(ip) becomes the second-order form, 
J(ip) ~ sin 2(p, where ip is the relative phase difference 1118 1. 
This is because a s-wave pairing is even under the mirror re¬ 
flection, while the A st ate is odd as shown in Table ED In 
addition, Yamakage et al. Ill 1811 pointed out that in the junction 


The search for Majorana fermion is currently one of the 
hottest issues in the physics community and various types of 
promising platforms have been theoretically proposed, as well 
as experimentally investigated during the last several years. In 
this review we focused on studies of superconducting doped 
topological materials that preserve the time-reversal-invariance. 
We began with a theoretical description of odd-parity pairing 
systems from a basic concept of symmetry invariances which 
are important to understand topological superconductors. The 
normal state of time-reversal-invariant topological supercon¬ 
ductors is characterized by the presence of Dirac fermions. We 
discussed the role of the Dirac fermions in the topological su¬ 
perconductors. To clarify the conditions for realizing topolog¬ 
ical superconductivity in doped topological materials, topolog¬ 
ical invariants and possible paring symmetry in the system are 
denoted. We continued with a summary of studies on the prop¬ 
erties of both normal and superconducting states in real ex¬ 
isting materials: Cu v Bi 2 Se 3 , Sni_ r In v Te, (Pbo.sSno.s) | Jn v Te, 
and Cu v (PbSe) 5 (Bi 2 Se 3 ) 6 . We found that the recent theoretical 
developments allow us to consistently interpret controversial re¬ 
sults aroused from different experimental methods. In this re¬ 
view, we shed light on the effects of disorder in the system. In 
particular, we discussed that the observed moderate suppression 
of T c of Cu v Bi 2 Se 3 with increasing the concentration of Cu, 
which behave as non-magnetic impurities in the superconduc¬ 
tor, can be explained by anomalous robustness of superconduc¬ 
tivity in the presence of nonmagnetic impurity scattering, which 
is peculiar to superconducting topological materials. Finally, 
in order to experimentally identify the existence of Majorana 
fermions in candidate platforms of topological superconduc¬ 
tors, we proposed future experiments to observe the anisotropic 
magnetic response of the candidate materials, as well as Joseph¬ 
son coupling effects. 
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